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SEA -LEVEL EVALUATION OF DIGITALLY IMPLEMENTED 
TURBOJET ENGINE CONTROL FUNCTIONS 
by Dale J. Arpasi, DavidS. Cwynar, and Robert E. Wallhagen 

Lewis Research Center 

SUMMARY 

The standard hydromechanical control system of a turbojet engine was replaced 
with a digital control system that implemented the same control laws. Three engine 
variables were controlled. The control required 1227 words of memory and 1. 4 milli- 
seconds for calculation. A detailed discussion of the digital control system in use with 
the engine is presented. The engine was operated in a sea- level test stand. 

The response under digital control for update intervals to 25 milliseconds is shown 
to be comparable to the response under hydromechanical control. At extended update 
intervals, a deterioration in transient response was noted. The use of digital compen- 
sation eliminated most of the response degradation. With compensation, closed- loop 
transient performance comparable to response under hydromechanical control was 
achieved at an update interval of 150 milliseconds. 


INTRODUCTION 

Digital control systems are receiving increased use in aircraft applications. Sys- 
tem complexities and interactions require the use of digital controls to obtain optimum 
performance under all operating conditions. An overall flight control system could in- 
clude controls for the airframe, the inlets, and the engines operating under the require- 
ments of mission control, navigation, and failure detection systems. Changes in re- 
quirements could change not only the performance criteria of the aircraft control sys- 
tems but also the control laws of these systems to arrive at optimum performance. The 
use of a digital computer in these flight systems permits maximum interaction and 
facilitates on-line modifications. 

The requirements of the digital computer system necessary to perform the aircraft 
control functions depend primarily upon the amount of calculation involved and the 



dynamics of the systems to be controlled. Being a discrete device, the computer sys- 
tem must be able to sample system measurements, to compute the control functions, 
and to put out commands at a rate fast enough so as not to impair the required closed- 
loop transient performance. For a sufficiently complex control law, the necessity of 
completing the prescribed calculations in the time interval specified by the update rate 
could lead to the need for very fast and costly computer systems. 

This report details the use of a digital computer to control a turbojet engine mount- 
ed in a sea- level static stand. The requirements for the digital system necessary to 
perform this task are indicated. A method is discussed for extending the control update 
interval to allow more computation time or to permit the use of a slower, less expen- 
sive computer to perform the same task. 

The digital system performs all functions of the standard hydromechanical control 
package, with the exception of afterburner control. Transient performance of the en- 
gine under digital control is presented for various update intervals. The update interval 
was extended and the deterioration in transient performance was noted. Digital com- 
pensation was incorporated into the measurements of the engine variables and resulted 
in improved transient performance. 


EXPERIMENTAL CONFIGURATION 

Digital engine control was carried out on a J85-13 engine operated in a sea- level 
static test stand. The standard hydromechanical control system of the engine was modi- 
fied to allow either electronic or hydromechanical computation of any control loop. 
Overrides, independent of the control, were incorporated to shut down the engine safely 
if engine overspeed or overtemperature limits are violated or if failures occur in the 
fuel, hydraulic, or electrical supplies. 


Engine Description 

The engine is illustrated in figure 1. Inlet air passes through inlet guide vanes into 
an eight-stage compressor. Air bleed doors are located on the third, fourth, and fifth 
compressor stages and along with the inlet guide vanes are referred to as the compres- 
sor variable geometry. The inlet guide vanes are mechanically coupled to the bleed 
doors, and the entire compressor variable geometry is driven by mechanically coupled 
hydraulic actuators located on opposite sides of the engine. The compressor variable 
geometry servomechanism control loop is scheduled to optimize compressor perform- 
ance over its complete range of operation. 
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I - Cutaway view oi J85 engine. 




A fuel control regulates fuel flow into the burner. The fuel flow passes through a 
pressurizing and drain valve and an overspeed governor before it enters the burner 
spray bars. The fuel and compressor variable geometry controls are located together 
in one of two hydromechanical packages located beneath the engine. 

Control of the variable exhaust nozzle area is provided by the other hydromechani- 
cal package. The output of this control is a mechanical position used to command three 
linear translating mechanical actuators coupled to the exhaust nozzle unison ring. 
Motion of the unison ring moves the nozzle leaves, thus varying nozzle throat area. 


Modification for Digital Control 

The standard hydromechanical engine control system was modified to permit selec- 
tive substitution of digital control for the standard hydromechanical controls of com- 
pressor variable geometry, fuel, and variable nozzle area and to permit easy switching 
between the two (for safety reasons). Hydromechanical computation of each control 
schedule was maintained, except for the compressor variable geometry. The hydraulic 
output from the standard hydromechanical control normally used to drive the compres- 
sor variable geometry actuators was inactivated. The compressor variable geometry 
schedule was implemented on a special-purpose analog computer. 

A block diagram of the experimental configuration is shown in figure 2. A servo- 
amplifier, of the type described in reference 1, was used to drive a servovalve to form 
the compressor variable geometry position loop. The servovalve powered the standard 
actuators; and their position, as measured by an engine- mounted linear potentiometer, 
was fed back to the servoamplifier to close the position loop. The command to the com- 
pressor variable geometry servoamplifier could be switched from either the analog 
computer schedule or the digital schedule. 

The main fuel control was left intact, with the exception of a four- way switching 
valve which was inserted in the engine fuel supply line just ahead of the overspeed 
governor. This allowed manual switching of engine fuel flow between the hydromechani- 
cal control and the digital control. The digital command scheduled the output of a fast- 
response fuel metering device (ref. 2) whose output was sent to the four- way switching 
valve. The selected fuel was ported by the valve to the engine while the unselected fuel 
was ported to drain through a dummy load. 

The standard exhaust nozzle control determines the position of a mechanical arm 
located on the exhaust nozzle control package. This arm is normally used to drive a 
power pack which in turn drives mechanical screw-jack actuators attached to the nozzle 
leaves through a unison ring. This system was modified by replacing the mechanical 
screw-jack actuators with hydraulic actuators and closing a servomechanism control 
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Figure 2. - Modifications to standard hydromechanical control system to permit switching to digital control system. 











loop around them in the same way as described for the compressor variable geometry 
control. The power pack was removed, and a linear displacement transducer was in- 
stalled to measure the position of the mechanical arm. The transducer output was 
scaled and linearized by an analog function generator so that the proper unison ring 
position command was presented to the exhaust nozzle servoamplifier. The command 
to the exhaust nozzle servoamplifier could be switched from either the analog function 
generator or the digital schedule. 

Signal conditioners were used to provide high-level signals (-10 V to +10 V) from 
the transducer measurements. Throttle position a was measured by attaching a 
rotating potentiometer to the shaft of the throttle input to the hydromechanical fuel con- 
trol. A static pressure tap attached to a high- response gage transducer was used to 
provide compressor discharge pressure Pg measurement. Engine speed N was 
measured by attaching a special gear to the engine power- takeoff and counting gear teeth 
with a magnetic pickup. The resultant frequency was converted to dc before being 
routed to the digital system. Thermocouples were used for measurement of both com- 
pressor inlet temperature Tg and turbine exit temperature Tg. An absolute pressure 
transducer was used to provide an ambient pressure Hq measurement for use in cor- 
recting Pg to absolute. 

The digital system ground was isolated from the experiment through isolation am- 
plifiers. This eliminated the possibility of ground loops and provided relatively noise- 
free measurements. 

CONTROL FUNCTIONS AND DIGITAL IMPLEMENTATION 

A block diagram of the digital system is given in figure 3 . It consists of a signal 
processing unit, an interface unit, a digital computer, and programming peripherals. 
The signal processing unit provides isolation of the signal inputs and calibration of the 
interface unit. The interface unit contains a high-level multiplexer, a single sample- 
and-hold amplifier, and a 13-bit digitizer to provide analog- to- digital conversion of the 
six measured signals. Digital- to- analog conversion of the computer commands is also 
accomplished in the interface through 13-bit digital- to- analog conversion units. The 
analog commands for compressor variable geometry position X v ^, fuel flow W^, and 
exhaust nozzle area Ag are obtained through separate digital- to- analog converters. 

The computer itself is a 16- bit machine with 16 384 words of magnetic core storage. 
The read-restore cycle time of the memory is 750 nanoseconds. Programming of the 
computer is done through a teletype and high-speed paper tape reader and punch. All 
control programming was done in assembly language to conserve core storage and ex- 
ecution time. 
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Figure 3. - Block diagram of digital control system. 


A complete description of this general-purpose digital system is given in refer- 
ence 3. Pertinent specifications of the signal processing unit, the interface units, and 
the computer are presented in table I. 


Compressor Variable Geometry Control 

The compressor variable geometry control schedules position as a function of 
corrected rotor speed N/^/6 and compressor inlet temperature Tg. In general, as 
WyfO increases, the compressor variable geometry is scheduled to close. Increasing 
Tg causes the variable geometry to open. The variable geometry includes the inlet 
guide vanes and the bleed doors. This control function was implemented on the digital 
computer by storing three curves of against N /^6 (each at constant Tg) in 

memory. Linear interpolation was used to compute X v ^ as a function of N /^6 and 
also to interpolate between lines of constant Tg. 

Fuel Control 

The fuel control commands a required ratio of fuel flow to static compressor dis- 
charge pressure W^/Pg as a function of throttle position a, engine speed N, and com- 
pressor inlet temperature Tg. The ratio is multiplied by the measured Pg, resulting 
in the required fuel flow. 
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TABLE I. - DIGITAL SYSTEM CAPABILITIES 


Digital computer 

Magnetic core memory size, words 

16 384 

Word length, bits plus parity 

16 

Memory cycle time, nsec 

750 

Add time, m sec 

1.5 

Subtract time, /isec 

1.5 

Multiply time, m sec 

4.5 

Divide time, /nsec 

8.25 

Load time, /usee 

1.5 

Store time, /usee 

1.5 

Indirect addressing 

Infinite 

Indexing 

Total memory 

Priority interrupts 

28 Separate levels 

Index registers: 


Independent 

1 

In conjunction with lower accumulator 

1 

Physical size, in. : 


Width 

24 

Height 

62 

Depth 

30 

Interval timers 

Complement 

2 

Accuracy, clock pulses 

±1 

Clock rates, kHz 

572, 286, 160, 143, 80, 71.5, 


40, 35.75, 20, 10 

Counter 

16- Bit binary 

Output 

Priority interrupt to computer 

Analog acquisition unit 

Overall sample rate (maximum), kHz 

20 

Resolution of digital data, bits 

12 (plus sign) 

Output code 

Two’s complement 

Number of channels 

64 

Input range, V full scale 

±10 

Input impedance, Mfi (shunted by 10 pF) 

10 

Maximum source resistance, 12 

1000 

Conversion time, Msec 

38 

Input setting time, Msec 

9 

Sam pie -and- hold aperture time, nsec 

500 

Safe input voltages, V 

±20 sustained 


±100 for less than 100 Msec 

Total error with calibration, percent 

0.073 




TABLE I. - Concluded. DIGITAL SYSTEM CAPABILITIES 


Frequency acquisition unit 

Number of channels 

10 

Nature of input 

Continuously varying or pulsatile 

Resolution of digital data, bits 

12 

Switch selectable clock rates, kHz 

20, 80, 100, 400, external 

Overall accuracy, bits 

±1 

Update rate 

Once per cycle of input frequency 

Maximum input frequency, kHz 

1 

Input amplitude range 

100 mV to 30 V peak to peak 

Analog output unit 

Total number of digital- to- analog 

26 

conversion channels 


Resolution (10 channels), bits 

12 (plus sign) 

Resolution (16 channels), bits 

11 (plus sign) 

Output voltage range, V full scale 

±10 

Output current (maximum), mA 

10 

Output impedance, Q 

<1 

Accuracy (12 bit), percent of full scale 

±0.1 

Accuracy (13 bit), percent of full scale 

±0. 05 

Slew rate, V/psec 

1 

Settling time for 10- V step to within 0. 05 

20 

percent of final value, psec 


Logical output unit 

Number of electronic switch outputs 

32 

Number of contact closure outputs 

32 

Maximum voltage, V 

30 

Maximum current, mA 

100 

Priority interrupt processor 

Number of channels 

10 

Input impedance, kfi 

47 

Input voltage range, V 

±10 

Comparator switching 

Trigger on rise or fall 

Comparator hysteresis 

Adjustable from 35 mV to 650 mV 

Comparator output, V 

+7 

Monostable multivibrator: 


Pulse width, psec 

. 3 

Pulse height, V 

+7 













A typical operating map of W^/Pg against N/y/0 is shown in figure 4. Droop 
lines (lines of constant a) are shown intersecting a normal steady- state operating line. 
A change in a, from to a^, at constant N/\/0 causes a step in W^/Pg. The 
W^/Pg command is limited by the acceleration schedules, which are functions of T2. 
As N increases, the W^/Pg command follows the proper acceleration schedule until 
it intersects the oig droop line, which it then follows to its steady- state value. The 
acceleration schedules are shaped to avoid turbine overtemperature and compressor 
stall during acceleration. 

In the droop- line calculation on the digital computer, desired steady- state values of 
Wj/Pg and N were computed as functions of a and Tg. Actual engine speed was 
subtracted from the desired speed to form a speed error. This error was multiplied 


by a proportional gain and added to the desired W^/P 
line output, (W^/Pg^ . The gain was proportional to 


g value to form the desired droop- 
N, giving the parabolic shape to 


the droop lines. 


In order to generate the acceleration schedules on the digital computer, eight 
curves of W f /Pg against N/-09 at constant Tg were stored in the computer memory. 
The measured Tg was compared to each of the eight values until the immediately lower 
value, T min , and the immediately higher value, T max , were found. Linear interpola- 
tion was used to find the W f /Pg value corresponding to N/y^0 at T min and T max< 
Linear interpolation was again used between these values to find the acceleration limit 



Corrected speed, N/V§, percent 

Figure 4. - Typical engine fuel control operating map showing acceleration and deceleration 
limits. 
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value, ( W f/ P 3 ) > corresponding to the measured T 2 . 

Upon completion of the droop and acceleration calculations, the smaller of 

and was selected. If the acceleration limit was chosen, a logic indicator V g 

was set as an acceleration indication to the exhaust nozzle control. A minimum limit 
on Wf/Pg was incorporated to provide a deceleration schedule. The result was then 
multiplied by the measured Pg to form the fuel- flow command. This command had a 
minimum limit at a value of fuel flow necessary to sustain combustion. For engine 
shutdown, a function was included which commanded zero fuel flow as the throttle was 
completely retarded. 


d 


Exhaust Nozzle Area Control 

The exhaust nozzle area Ag was controlled primarily as a function of a. Certain 

engine conditions caused an override of the throttle- commanded area. If the engine was 

operating on the fuel control's acceleration limit (V set), the minimum area was made 

s 

equal to a fixed area, Ag a , to allow faster acceleration. At low compressor discharge 
pressure the area was also limited from closing to less than a function of Pg , Agp, to 
provide compressor stall protection. 

The exhaust nozzle control also provides turbine overtemperature protection by 
overriding the normal throttle- scheduled area with an integral- proportional regulation 
of turbine exhaust gas temperature Tg. This regulation was implemented on the digital 
computer by taking the Z- transform of the control function used in the standard hydro- 
mechanical control. The resultant Tg override area Agj. was limited to its full- 
closed value whenever Tg was less than its reference value. A maximum selector was 
used to pick the larger of the normal area command or Ag t> Consequently, if Tg ex- 
ceeded its reference value, the nozzle would open and reduce Tg. The reference value 
was made a function of a nominal Tg limit, Tg, and N. 

Input-Output and Timing 

The three controls (compressor variable geometry, fuel, and exhaust nozzle), im- 
plemented as described, require the measurement of six variables - a, N, Tg, Tg, Pg, 
and Hg. A high-level (-10 V to +10 V) multiplexer, a sample- and- hold, and a digitizer 
under block data transfer control (ref. 3) were used to take sampled values of these 
measurements directly into computer memory on a cycle- stealing basis. That is, upon 
computer command, channel transfer to the multiplexer and data transfer from the 
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digitizer were accomplished automatically without interrupting the computer for more 
than one memory cycle per word transferred. The sample- hold- digitizing processes 
operated at a rate of 20 000 words per second (50 /rsec/word). Acquisition of the six 
measured variables, therefore, required 300 microseconds. A priority interrupt was 
issued by the block transfer control to the computer when all six words were in memory. 

A 13- bit digital- to- analog converter was assigned to each of the control outputs 
(X V g, Wj, and Ag). Upon completion of the control computation, these commands were 
issued from the computer and required 6. 75 microseconds per word transferred. 

A timing diagram giving the event sequence over an update interval AT is presented 
in figure 5. The control process was sequenced through the use of an interval timer. 


Timer 


Sample- 

complete 


Time, t, 
psec: 


Reinitialize timer 

Idle during 

Scale sampled parameters; 

Idle until next 

* ano sxan measure • 
ment sampling 

sampling 

“""compute control and output 

timer interrupt - * 


Timer 

interrupt 


7.5 


307.5 


-1408 


AT 


Depends on 
computational 
route (maximum 
time indicated) 


Figure 5. - Timing diagram for typical update interval. 


This timer accepts a digital word from the computer and counts it down at a fixed rate. 
When the countdown reaches zero, the timer issues a priority interrupt to the computer, 
indicating that a fixed amount of time has elapsed. The magnitude of the word dictates 
the length of the update interval AT. Upon receiving the timer interrupt, the computer 
reinitializes the timer for the next update interval and begins sampling the six measured 
variables and taking the digitized values into assigned memory locations. This requires 
about 7. 5 microseconds. The computer then idles for 300 microseconds during the 
digitizing period. During idle periods, the computer is available for processing other 
routines if so required. 

When the measurements are in memory, the priority interrupt from the block data 
transfer control occurs and the control computation begins. This computation time de- 
pends on the route taken through the program and averages about 1. 1 milliseconds. Upon 
completion of the computation the manipulated variables are updated and the computer 
idles until the next priority interrupt from the interval timer. 

Calculation times and memory requirements for each subroutine and the complete 
control calculation are presented in table n. This table provides time information in 
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TABLE II. - CONTROL COMPUTATION MEMORY AND TIMING REQUIREMENTS 
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terms of memory cycles (to allow the reader to estimate calculation time for other com- 
puters) and in terms of milliseconds based on the computer used in this experiment. The 
complete engine control program requires 1227 words of memory and a maximum of 
1. 4 milliseconds of time including sampling. The maximum calculation time is based on 
manufacturer published times for the computer instruction set. 


RESULTS AND DISCUSSION 

The control update interval AT determines how often the engine variables are sam- 
pled and how often the commands are updated. The minimum possible AT is set by the 
sum of the sampling time and calculation time of the digital system. Table n indicates 
that for this system the maximum calculation time is 1. 15 milliseconds and the time nec- 
essary to sample the six measured variables is 0. 3 millisecond. The minimum AT 
possible with this control is, therefore, 1. 5 milliseconds. 

Engine transient response dictates the maximum possible update interval. Since 
the engine is a nonlinear device, it would not suffice to evaluate the effects of increased 
AT by using linear stability analysis at a single operating point. For this reason, 
throttle position steps from idle (approx. 50-percent N) to military (100- percent N) were 
used to provide the engine transients for update interval evaluation purposes. 

The dashed lines in figure 6 illustrate the engine response under the modified hydro- 



Figure 6. - Comparison of digital and hydromechanical controls for a throttle step from idle 
to military. 
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mechanical control (fig. 2) for a throttle slam from idle to military. Hereinafter, the 
modified hydromechanical control of figure 2 will be referred to simply as the hydro- 
mechanical control. Rotor speed N, turbine discharge temperature Tg, compressor 
exit static pressure Pg, and gross thrust are plotted against time. The same 
figure gives the response of the engine under digital control (solid lines) for an update 
interval AT of 2 milliseconds. The digital response does not exactly match the modi- 
fied hydromechanical response. It exhibits a somewhat shorter rise time and less 
overshoot in rotor speed, while allowing more overshoot in turbine exhaust gas temper- 
ature. These differences were traced to slight schedule differences in the main- fuel 
and variable exhaust nozzle control systems. The digital schedules were taken from 
and checked against information received from the manufacturer. 

Figure 7 is a plot of compressor pressure ratio against corrected speed during the 
same throttle slam from idle to military. The response under hydromechanical control 
is approximately equivalent to the response under digital control. The apparent dis- 
crepancy at the high end of the transient is due to differences in the exhaust nozzle con- 
trol. Under hydromechanical control the nozzle closes earlier, causing a faster rise in 
pressure ratio and increased overshoot in corrected speed. 

The response of the engine to the large-scale transients may be evaluated in terms 



Figure 7. - Effect of digital and hydromechanical controls on 
compressor pressure ratio for a throttle step from idle to 
military. 
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of the time required to move from one operating level to another and how well it re- 
mained within the prescribed operating limits of the engine. These limits are engine 
overspeed, turbine overtemperature, compressor stall, and blade vibrations. The 
engine was not instrumented to measure blade vibrations. However, engine vibration 
was monitored by panel meters and was well within prescribed limits for all conditions 
discussed in this report. The engine response under the modified hydromechanical con- 
trol was used to provide a conservative estimate of the operational limits. 

Comparisons of the responses in figures 6 and 7 illustrate that under digital control 
at the 2-millisecond update interval the engine response equals or surpasses the re- 
sponse under hydromechanical control. The possible exception is the turbine exhaust 
gas temperature overshoot. This overshoot is a function of the integral and proportion- 
al gains in the Tg override of the variable exhaust nozzle control; it could have been 
trimmed to exactly match the hydromechanical characteristics. However, since the 
real temperature limit is that of the turbine blades, this slight overshoot in gas tem- 
perature was considered negligible . 


Effects of Control Update Interval 

The effects of update interval on engine response are shown in figures 8 and 9, for 



Figure 8. - Effects of update interval on engine response for a throttle step from idle to military. 
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Figure 9. - Effect of update interval on compressor pressure 
ratio for a throttle step from idle to military. 


update intervals of 2, 25, 50, 100, 150, and 200 milliseconds. For simplification pur- 
poses, only gross thrust and rotor speed are illustrated in figure 8 and only the 2-, 

25-, and 2 00- millisecond intervals are illustrated in figure 9. The engine response 
was stable throughout this range. However, as the update interval was increased, the 
engine response deteriorated in terms of rise time; and between 25 and 50 milliseconds 
the overshoot in rotor speed exceeded the 2-millisecond reference response. Figure 9 
indicates that the pressure ratio transient response at a AT of 200 milliseconds did 
not appreciably exceed the 2- millisecond reference response. 

From the data of figures 8 and 9, it can be stated that the digital control with a AT 
of 25 milliseconds was indistinguishable from the control with a AT of 2 milliseconds. 
At extended intervals the response deteriorated (although not extremely), and violation 
of the engine operational limits was possible. 


Use of Digital Compensation to Extend Update Interval 


The deterioration in response with increased update interval is due to the fact that 



the sampled values of the measurements at the start of an interval are an increasingly 
poor approximation to the values of the variables over the interval. The value of a 
variable in the n** 1 time interval e n (t) can be approximated by a power series (ref. 4) 
of the form 


e„(t) = e(n AT) + e(n AT)(t - n AT) + i- e(n AT)(t - n AT) 2 + . 
n 2! 


( 1 ) 


where e(n AT) is the sampled value of the variable at time n AT, and e and e are 
the first and second derivatives, respectively. (All symbols are defined in the appen- 
dix. ) Time t goes from n AT to (n + 1)AT. If AT is appreciable and the measured 
variable is changing rapidly, much information on the state of the variable is missing 
from the sampled value. 

The deterioration in response with increased AT can be further explained by 
studying the frequency response of a digital sampler and zero- order hold shown in fig- 
ure 10. The frequency axis is shown as a function of the sampling frequency f = 1/(AT). 
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Figure 10. - Normalized frequency response of a sample- 
and-hold. G(sl HI -e' ATs )/at s. 


As long as f g is large compared to the engine frequencies f, the sample- and- hold adds 
little to the closed- loop dynamics (amplitude so. 98, phase s- 17. 5°). Decreasing f 

s 

(increasing f/f g beyond 0. 1), however, introduces appreciable attenuation and phase 
shift. If we could compensate for this attenuation and phase shift, we could extend the 
update interval (decrease f g ) while maintaining good engine response. 

A typical sample- hold- compensation operation is depicted in the block diagram of 
figure 11. The variable e m is any engine variable sampled and held to produce the 
discrete variable e^. Digital compensation of e h results in the discrete variable e c - 
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Figure 11. - Block diagram of digitally compensated sample-and-hold of a typical engine variable. 


In an effort to provide some phase lead, the compensator was made to be a linear func- 
tion of the digital estimations to the first and second derivatives of e h and is given by 

e c (n AT) = e h (n AT) + Kj AT e h (n AT) + Kg AT 2 e h (n AT) (2) 

The derivative gains and K 2 were variable and adjusted experimentally to improve 
the engine response to throttle steps. 

The digital estimations of derivatives are 

6 >,< n aT > - ± aT > - e h 0" - d at ]} < 3 > 

e h (n AT) =-J— [e h (n AT) - 2e h [(n - 1)AT] + e h [(n - 2)AT]} (4) 

(AT) 2 

Substituting equations (3) and (4) into equation (2) produces the digital equation 

e c (n AT) = (1 + Kj + K 2 )e h (n AT) - (Kj + 2K 2 )e h [(n - 1)AT] + K 2 e h [(n - 2 )At] (5) 

In order to observe the effects of compensation in the frequency domain, equation (5) 
may be rewritten in Laplace transform notation as 

= (1 + K 1 + K,) - (Ki + 2K ? )e" ATs + K 9 e" 2ATs (6) 

e h (s) 12 12 *2 

The net result of including this digital compensation is given by the transfer function 
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[(1 + K 1 + K 2 ) - (K t + 2K 2 )e“ ATs + j^^ATsJ 


( 7 ) 


e c (s > 

e m ( s ) 


1 - e 


-ATs 


AT s 


A frequency domain plot of this transfer function is shown in figure 12 for = 1 and 
K 2 = 0. 5. When this response is compared to that of the uncompensated sample- and- 
hold shown in figure 10, the large amount of phase lead added by the compensation be- 
comes apparent. However, it should also be noted that substantial gain is also inserted 




Figure 12. - Normalized frequency response of a compen- 
sated sample-and-hold. GIs) ■ [il -e' AT s )/AT s] 
<2.5-2e' ATs + 0.5e' ZATs I. 


into the control loop by the compensation, which tends to reduce the stability margin of 
the closed loop. 

To demonstrate the effects of this type of compensation on the engine at extended 
update intervals, the relation given by equation (5) was programmed into the measure- 
ment of N, Pg, and Tg. That is, each of these variables was sampled and held as 
before, but before being acted upon by the control program they were each compensated 
according to equation (5). The amount of compensation (K^ and K 2 ) was increased 
from zero at a AT of 150 milliseconds and the response of the engine to throttle slams 
was observed. 

An example of the effects of compensation on the engine response to throttle slams 
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is given in figures 13 and 14. The sampled values of N, Pg, and Tg were compen- 
sated according to the equations 


N c (n AT) = 2 . 5N h (n AT) - 2N h [ 

(n - 1)AT] + 0. 5N h [(n - 2)AT] 

(8) 

P 3c< naT >= 2 - 5P 3h< naT >- 2P 3h 

[(n - 1)AT] + 0. 5P 3h [(n - 2)AT] 

(9) 

T 5c (nAT)=2.5T 5h (nAT)- 2T 5h 

[(n - 1)AT] + 0. 5T 5h [(n - 2)AT] 

(10) 


For comparison purposes the engine responses without compensation at update intervals 
of 2 milliseconds and 150 milliseconds are also shown in this figure. The use of com- 
pensation essentially allowed the response at 150 milliseconds to match the uncompen- 
sated response at 2 milliseconds. It was found, however, that any further increases in 
the derivative gains and Kg resulted in instability at high- rotor- speed operation. 
This instability was manifested in terms of diverging oscillations in fuel flow. It was 
also determined that the response with this compensation was very sensitive to changes 
in compressor inlet temperature Tg (Tg for figs. 13 and 14 was approximately 311 K 
(100° F) at military operation). Unfortunately, documentation of the effects of tempera- 
ture was not possible because of the inability to regulate this variable in the test facility. 

It should be remembered that responses for sampling intervals up to 25 millisec- 
onds were essentially equivalent to that at 2 milliseconds. Therefore, operation at 
150 milliseconds as opposed to 2 milliseconds results in an increase in computer effi- 
ciency of 6 to 1. At 150 milliseconds only 1/1 00th of the interval is being used in the 
control computation. The remainder of the time may be used to provide other aircraft 
control and monitoring services. In another regard, the use of compensation allows the 
use of slower and less expensive computing hardware while maintaining the response of 
a fast machine. The compensation described is by no means optimum or unique but in- 
dicates the advantages of the use of digital compensation in aircraft control applications. 


CONCLUSIONS 

A digital control system to replace the hydromechanical control of compressor 
variable geometry, main fuel flow, and variable exhaust nozzle area of a J85-13 turbo- 
jet engine requires a computer with slightly less than 1500 words of core storage. The 
entire digital control process can be calculated in approximately 1500 computer memory 
cycles. The repetitive update interval may be extended to 25 milliseconds without ap- 
preciable degradation in closed- loop response. 
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Digital compensation, making use of weighted approximations to the derivatives to 
incorporate prediction into the sampled engine variables, allows extension of the maxi- 
mum update interval to 150 milliseconds with no appreciable loss in engine response. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 12, 1972, 

764-74. 
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APPENDIX -SYMBOLS 


2 

exhaust nozzle area, cm 

o 

exhaust nozzle area during acceleration, cm 

exhaust nozzle area as function of pressure, cm 

exhaust nozzle area as function of temperature, 

any engine variable 

first time derivative of e 

second time derivative of e 

gross thrust, N 

frequency, Hz 

sampling frequency, Hz 

transfer function (appropriate dimensions) 

2 

ambient total pressure, N/cm 

2 

compressor inlet total pressure, N/cm 

o 

compressor exit total pressure, N/cm 
first derivative gain, dimensionless 
second derivative gain, dimensionless 
percent rotor speed 
percent corrected rotor speed 
integer 

2 

compressor discharge static pressure, N/cm 
Laplace variable, sec - * 
temperature, K 

compressor inlet temperature, K 

turbine exhaust temperature, K 

sampling interval, sec 

time, sec 

logic indicator 

fuel flow, kg/hr 



X variable geometry position, percent 

v o 

o; throttle position, deg 

e exponential function 

Subscripts: 
a acceleration 

c digital compensator output 

d droop 

h sample -and- hold output 

m measured 

max maximum 

mil military 

min minimum 
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